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A new method for predicting inflationary cosmological perturbations, based on the Wentzel- 
Kramers-Brillouin (WKB) approximation, is presented. A general expression for the WKB scalar 
and tensor power spectra is derived. The main advantage of the new scheme of approximation is that 
it is valid even if the slow-roll conditions are violated. The method is applied to power-law inflation, 
which allows a comparison with an exact result. It is demonstrated that the WKB approximation 
predicts the spectral indices exactly and the amplitude with an error lower than 10%, even in regimes 
far from scale-invariance. The new method of approximation is also applied to a situation where 
the slow-roll conditions hold. It is shown that the result obtained bears close resemblance with the 
standard slow-roll calculation. Finally, some possible improvements are briefly mentioned. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



Observations of the inhomogeneities and anisotropies 
of the Universe's large scale structure and of the cos- 
mic microwave background (CMB) radiation are one of 
the keys to decipher the cosmic evolution. According to 
the favoured cosmological scenario, fluctuations of mat- 
ter and space-time are generated during an early epoch 
of inflation. The ability to predict the power spectra of 
cosmological perturbations for a large variety of differ- 
ent inflationary models is of the utmost importance. So 
far, all predictions that come up to the level of accuracy 
needed for CMB experiments such as Archeops, MAP or 
Planck, have to rely on the slow-roll approximation, sub- 
ject to an infinite number of convergence conditions, or 
on numerical integration Jj. In this paper we present a 
new method that overcomes the mentioned restrictions 
(but has its own, very different ones). Our method is 
applicable to models that cannot be described by the 
slow-roll approximation (e.g. the 'problem of a large 
slow-roll parameter 77' of some supersymmetry-inspired 
inflationary models) and it can be applied to small scales 
that leave the horizon at times close to the end of infla- 
tion, where the slow-roll approximation necessarily has 
to break down. 

The study of scalar and tensor fluctuations can be re- 
duced to the study of a single variable during inflation, 
usually denoted by /J,g or /z^ [2| • The corresponding equa- 
tions of motion are similar to the Schrodinger equation 
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for a one-dimensional system: 

+ [k' - uiv)]^^ = + = , (1) 

where, in the cosmological context, the conformal time 77 
plays the role of the radial variable r in atomic physics. 
In general, this equation cannot be solved exactly un- 
less very simple potentials U{ri) are considered and one 
has to rely on a method of approximation. However, we 
have inherited a very powerful method to deal with more 
complicated potentials from the pioneers of quantum me- 
chanics, namely the Wcntzcl-Kramers-Brillouin (WKB) 
method of approximation. 

In this article, we apply this method to the equation 
governing the evolution of cosmological perturbations. It 
is quite astonishing that for more than twenty years since 
the discovery of cosmological inflation, to our knowl- 
edge, nobody has been able to apply successfully the 
WKB approximation to the prediction of cosmological 
perturbations. This is probably related to the fact that 
the growth of the perturbations on superhorizon scales, 
which is equivalent, in the quantum- mechanical version 
of the theory, to particles creation, is often described in 
terms of the breakdown of the WKB approximation. In 
this article, we demonstrate that this is true only if the 
WKB approximation is used naively. In fact, the problem 
bears a close resemblance with the situation discussed by 
atomic physicists at the time quantum mechanics was 
born. We therefore make a historical digression to illus- 
trate the main point of the present article. 

The subject debated by the atomic physicists at the 
beginning of the twentieth century was the application 
of the WKB approximation to the motion in a central 
field of force and, more specifically, how the Balmer for- 
mula, for the energy levels of hydrogenic atoms, can be 
recovered within the WKB approximation. The effective 



2 



frequency for hydrogenic atoms is given by (obviously, 
in the atomic physics context, the wave equation is not 
a differential equation with respect to time but to the 
radial coordinate r) 

-"W = ^(i^+— j- (2) 

where Ze is the (attractive) central charge and £ the 
quantum number of angular momentum. The symbol E 
denotes the energy of the particle and is negative in the 
case of a bound state. The calculation of the energy lev- 
els was first addressed by Kramers 3] and by Young and 
Uhlenbeck Q . The method employed was Sommerfeld's 
quantization rule 

^tj(r)dr = 27r^7i+ , (3) 

where the integral must be calculated along the classical 
trajectory given here by an ellipse. The result obtained 
by those authors reads 

In addition, the WKB wave function for small values of 
r was found to behave as ~ rl^^^+^^l ^ +i/2_ xhis last re- 
sult follows straightforwardly if one retains only the last 
term in Eq. (^J, i.e. uj{r) ~ i[i{(+l)]^^'^/r, as required in 
the limit r ^ 0. It was noticed by Young and Uhlenbeck 
that the two previous results are not satisfactory. Indeed, 
experimental results indicate that the factor i{£ + 1) in 
Eq. should be replaced by {£ + 1/2)^: "It is apparent 
that we do not obtain the familiar Balmer formula ex- 
cept in the limit of large £ values. If, however, we replace 
£{£+!) by (^-t-1/2)^ the correct result is obtained^ A sim- 
ilar remark is made in Ref. 4] concerning the behaviour 
of the wave function near the origin: "The approximate 
solution based upon the negative exponential behaves for 
r very small like ' +1/2^ Again the replacement 

of £{£ + 1) by {£ + 1/2)^ leads to an agreement between 
the approximate and exact solutions." The solution advo- 
cated by Young and Uhlenbeck was therefore to replace, 
by hand, ^(£-1-1) by (£+1/2)'^ without further theoretical 
justifications. 

In 1937 the problem was considered again by 
Langer In the section devoted to the radial wave 
equation, he notices that the WKB approximation was 
applied to the Schrodinger equation Q by Kramers, 
Young and Uhlenbeck only because of its resemblance 
with the equation of an oscillator: "Insofar as the writer 
is aware, the degree of this resemblance has hitherto been 
regarded without exception as sufficient for assuming the 
out-and-out applicability to the equation (17) of the for- 
mulas deduced in the foregoing discussion, especially of 
the formula (14)," where Eqs. (17) and (14) in Langer's 
paper correspond to Eq. and the WKB connection 
formula. He also remarks that no theoretical explanation 



has been given to explain why the wrong factor £{£ -I- 1) 
is obtained: "This "failure" of the WKB method, i.e. to 
the extent that the change in question is requisite, has 
been generally verified both in studies of attractive and 
repulsive fields. No explanation of it seems to have been 
given, though that can be done very simply to the follow- 
ing effect. The fault lies not in the method but in the 
application of it." Then, Langer shows that the WKB 
approximation breaks down, at small r, for an effective 
frequency given by Eq. |(2Jl and, in addition, he suggests 
a method to circumvent this difficulty The method con- 
sists in changing the variables according to 

r = e^, /i = e^/2u. (5) 

The resulting equation for u has again the form of Eq. 
but the effective frequency is now given by 

c.^(x) = |j(i?e^- + ZeV)-(^+0'. (6) 

As demonstrated by Langer, the crucial point is 
that, now, the WKB approximation, applied to the 
Schrodinger equation for u{x), no longer breaks down. 
Therefore, it is now possible to utilize the WKB formal- 
ism and to use the usual WKB approximation. Trans- 
forming back to the original variables the correct factor 
(^-fl/2)^ appears, instead of the wrong term [£(^-1-1)]^/^. 
This allows Langer to conclude that his method "... is 
seen to amount formally to precisely the replacement of 
£{£+1) by (£+1/2)2". 

Let us now come back to the theory of cosmological 
perturbations. The analogy with the previous situation 
is that large values of r correspond to subhorizon scales 
(A: 1 77 1 S> 1) whereas small values of r correspond to super- 
horizon scales {k\ri\ <C 1). On subhorizon scales the mode 
function oscillates (w ~ k) and the WKB approximation 
can be applied in this regime without any difficulty. On 
superhorizon scales the perturbations do not oscillate, 
but this is not a problem in itself. The WKB approxima- 
tion can be used in a situation where the potential energy 
dominates and leads to exponentially decaying/growing 
solutions (e.g. in the WKB treatment of the a decay of 
a radioactive nucleus). Nevertheless, difficulties arise if 
one tries to find the WKB approximation to the mode 
equations of cosmological perturbations. 

On superhorizon scales a problem arises in the naive 
WKB approach a la Kramers, Uhlenbeck and Young from 
the form of the effective potential U{ri). It turns out 
to be the very same problem as encountered by atomic 
physicists. The resemblance is most explicit in the par- 
ticular case of power-law inflation where the potential 
U{ri) (X 77"^ is exactly of the form given in Eq. ||2J). The 
usual conclusion that the WKB approximation is vio- 
lated for superhorizon scales is in fact not true, provided 
the cosmological version of the transformation jS)) can 
be found. This opens the possibility to employ the WKB 
approximation to calculate the spectrum of inflationary 
cosmological perturbations in the superhorizon limit. In 
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the cosmological context, the WKB approximation is a 
new method for computing the power spectra, different 
from the more traditional slow-roU approximation. 

This article is organized as follows. In the next sec- 
tion, we briefly recall the basic results of the inflationary 
cosmological perturbations theory. In the third section, 
we show how the WKB approximation can be applied 
to inflationary cosmological perturbations in the general 
case, i.e. for any inflaton potential. The only restric- 
tion is that we focus on single inflaton field models. We 
demonstrate that the calculation of the spectrum reduces 
to the calculation of a single quadrature. In the fourth 
section we apply the WKB approximation to power-law 
inflation. This model is exactly solvable so that we can 
compare the exact result with the WKB prediction and 
also with the standard calculation based on the slow-roll 
approximation. This allows us to study the precision of 
the WKB approximation in the cosmological context. In 
particular, we demonstrate that, in the regime where the 
slow-roll approximation breaks down (far from scale in- 
variance), the WKB approximation is still valid. This 
opens the possibility to calculate the power spectra of 
more realistic models violating the slow-roll conditions. 
In the fifth section, we apply the WKB approximation to 
the general (first-order) slow-roll situation. We compare 
the result with the standard slow-roll calculation for the 
amplitudes and the spectral indices. We show that the 
spectral indices are very well predicted while the error 
in the amplitude is of the order 10%, i.e. less accurate 
than the prediction of the slow-roll approximation for 
this quantity. Finally, in the sixth section, we give our 
conclusions and indicate what the possible extensions to 
the present work are. 



equation of motion for the mode functions fi should be 
solved with the following initial condition 
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(8) 



where r]i is an arbitrary time at the beginning of infla- 
tion and rripj is the Planck mass. This initial condi- 
tion corresponds to the fact that, initially, the modes 
are subhorizon and, therefore, do not feel the curvature 
of space-time. As a consequence, they are described by 
plane waves. If the initial quantum state is the vacuum 
state, then the statistical properties of the perturbations 
are entirely characterized by the two-point correlation 
function, i.e. by the power spectrum. The dimension- 
less power spectra of scalar and tensor fluctuations are 
calculated to read 



87r2 



Vh = 



(9) 



The spectral indices and their running are defined by 
the coefficients of Taylor expansions of the power spectra 
with respect to In fc, evaluated at an arbitrary pivot scale 
/c*. To be more specific: 



1 



dlnVc; 



dlnfc 
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k=k. 
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fe=fc. 



are the spectral indices and the two following expressions 



d2 In -Pc 
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(dlnfc)2 



d^ InVh 
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(dlnfc) 2 



(11) 



k=k. 



II. GENERAL EXPRESSIONS FOR 
INFLATIONARY PERTURBATIONS 



define the "running" of these indices. In principle, we 
could also define the running of the running and so on. 



A. Power spectra 



B. Inflationary parameters 



The evolution of the cosmological perturbations in 
a single field inflationary scenario is governed by the 
Schrodinger-like equation The effective time- 

dependent frequency is given by the general expression 



z 



(7) 



where the primes denote a derivatives with respect to 
conformal time and where — a\j2 — aa" /a''^ for scalar 
(density) perturbations and = a for tensor perturba- 
tions; a(?7) is the scale factor of the Friedmann-Lemaitre- 
Robertson- Walker metric. For density perturbations, 
the quantity /i is given by /ig ~ ^SzgC, where C, is 
Bardeen's hypersurface- independent quantity Q, which 
for adiabatic perturbations is conserved on superhorizon 
scales. For tensor fluctuations, one simply has /i.^, = z^h, 
where h is the amplitude of gravitational waves. The 



In order to investigate the predictions of cosmological 
inflation in a model-independent manner, it is useful to 
define a set of functions that uniquely characterize the 
state of the Universe at any given moment of time. It 
is usefuU to describe the evolution of the Universe by 
a set of flow equations 0, In this article we make 
use of the set of horizon flow functions, introduced in 
Ref. 0. The zeroth horizon flow function is defined by 
eo = H{Ni)/H{N), where H = a/a is the Hubble rate, a 
dot meaning derivative with respect to cosmic time. In 
this expression, N is the number of e-folds, N = ln(a/ai), 
after an arbitrary initial time. The hierarchy of horizon 
flow functions is then deflned according to 
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dA^ 



71 > 0. 



(12) 



Inflation takes place for ei < 1. It is interesting to es- 
tablish the link between the horizon flow functions and 
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the so-called slow-roll parameters, the first three being 
defined by 



^ 

Hip 
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(13) 
(14) 
(15) 



where ip is the inflaton field and V{ip) is the inflaton po- 
tential. The first two slow-roll parameters must be small 
during inflation, as can be seen from their definition. The 
parameter e is just the inflaton kinetic energy to inflaton 
total energy ratio, whereas 5 is the acceleration to speed 
ratio. The set {e, S, £,} is linked to the horizon flow func- 
tions {e„} by 



1 



(16) 



The set {e, S, ^} or equivalently the set {ei, £2, £3} plays 
a special role, since the effective potential U{ri) = z" j z of 
Eq. can be expressed exactly in terms of these three 
functions only. One finds 



= a^H^{2~ei). 



1 ,2 



^e2e3),(17) 
(18) 



Let us emphasize that these expressions are exact, no 
approximation has been used so far. 



III. WKB POWER SPECTRA 

We now turn to the derivation of the WKB inflationary 
power spectra for scalar and tensor perturbations. We 
define the quantity /i^^B by the expression 



MwKB 
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tjl/2(fc^,y) 



,±i/''w(fe,T)dT 



(19) 



The mode function /i^^e represents the leading order 
term of the semiclassical expansion. It satisfies the fol- 
lowing differential equation 



u;\k,r^)-Qik,rj) 



MwKB(fc,»/) = 0, (20) 



where the quantity Q{k,r]) is given by 



CO 

2lo 



(21) 



Therefore, the mode function /^wkb(^i'7) given in 
Eq. H19|) is a good approximation of the actual mode 
function fi(k,r]), if the following condition is satisfied 



Q 



< 1 



(22) 



On subhorizon scales, lo k, which implies Q — 
and therefore the condition is satisfied. On supcrhorizon 
scales, it is not possible to show that the previous condi- 
tion is always violated, since it depends on the shape of 
the effective potential U{ri). However, in the cosmological 
context, the problem discussed in the introduction gen- 
erally arises. We explicitly show this in the next section 
for the case of power-law inflation. Here, we demonstrate 
the problem for the more generic case of slow-roll infla- 
tion. Using the expression for the effective frequency in 
terms of the horizon flow functions, one has 



Q 
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~ 8 

T 



0(ei,e2). 



(23) 



The corrections, linear in the horizon flow functions, 
could easily multiply the term 1/8 = 0.125 by a factor 2. 
Therefore, the number \Q/uj^\ is indeed not greater than 
1 but, at the same time, not very small with respect 
to 1. Thus, the WKB approximation is not a good ap- 
proximation in this regime. In order to have an effective 
potential that permits to use the WKB approximation, 
we introduce the following transformation 



^Inf^), .^(l-eO^/V/V 



(24) 



This is the cosmological counterpart of Langer's transfor- 
mation given in the introduction; see Eq. |SJ). However, 
let us remark that there exist differences, in particular 
the factor (1 — ei)^/^ in the definition of the function u. 
It is also interesting to notice that this transformation 
has already been considered in the literature in Ref. 
for different reasons. Then, for density perturbations, 
the equation of motion takes the form 



d^^s(a:) 



1/3- 
4Vl 



ei 



ei 



(3 - 2£i)e2 _ (1 - 2ei)e2e3 _ (1 - 4£i)e^ 
2(l-ei)2 2(l-ei)3 4(1 - ei)4 



Itg = 0, 



(25) 



whereas, for gravitational waves, we obtain 



d^u^ {x) 
dx^ 



-(- 

4l 1 



ei 



ei 



eie2 



eie2e3 (2 -I- ei)eief 



2(l-ei)2 2(l-ei)3 4(1 -ei)^ 



u,j.(x) ~ 0. 



(26) 
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In the following lo^ and u J^ de note the expressions in the 
square brackets of Eqs. and (j^ . The WKB ap- 

proximation can now be applied to these equations. In 
order to demonstrate this, let us recalculate, on super- 
horizon scales, the quantity \Q/ijj^\ for the two previous 
expressions. One finds for density perturbations 



Q 








S 


27 \ 



6162- 



+Oiet), (27) 



whereas for gravitational waves, one has 



Q 



27 



(eiel + £16263) + 



(28) 



It is obvious that the quantity jQ/w^j is now very small 
for any inflaton potential satisfying the slow-roll condi- 
tions, since this is a third-order quantity. As a conse- 
quence, the WKB approximation is now also valid on 
superhorizon scales. 

Let us apply the WKB approximation to Eqs. H25|) 
and H26|l . Let be the value of the variable x such that 
ti;^(x*) — 0, i.e. is the classical turning point. We 
define region I as the region such that uj'^{x) > and 
region II the region where uj'^{x) < 0. Let Xi be the 
initial "time" at which the normalization is performed 
(i.e. corresponding to 771). At some point x, deep in 
region I, the WKB approximation of the function u{x) 
can be written as 



u,{x) 



A 



;i/2(a;) 



exp 



±i 



uj{y)dy 



(29) 



where A is a constant of normalization. Deep in region I, 
we have uj{y) ~ e~^/ (1 — ei) and straightforward manip- 



ulations show that the corresponding fi^ (rj) behaves as in 
Eq. IHl , provided that one chooses the minus sign in Eq. 
(|^ and that A is given by 



A 



(30) 



In region II, the effective frequency is pure imaginary 
and can be written as uj(x) = i|a;(x)|. Let X{ be the final 
point at which we evaluate the solution deep in region II. 
Then, at any point x in region II, the WKB solution is 
given by 



u,,{x) 



\co{x)\ 



-exp 



\u;{x)\ 



-exp 



k(2/)|dy 



(31) 



where C± are constant. The goal is now to connect these 
two constants with the constant A. At the turning point 
X = Xi,, the WKB approximation breaks down and one 
has to use the usual WKB procedure and approximate 
the potential in this region by a straight line such that 

u}^{x)~-a{x~x^), (32) 
where a = —[duj'^{x)/dx]{x — x^,) > 0. The solutions of 
Eqs. I|25|l and (|26|l are given in terms of Airy functions 
of first and second kinds 



u{x) = BiAi{s) + B2Bi{s), 



(33) 



where s = a^^^{x — x„) 



The next move is to use the asymptotic behaviour of the Airy functions to calculate the relation between Bi and 
B2 with A on one hand and C± on the other. In region I, for a value of x not too far from x*, Eq. H33|) can be written 
as 



u,{x) 



a 12 



\X -3^ * 



i2 1 



(i?2 - iBi) exp — a2 - x^\^ + — ] + {B2 + iSi) exp -—a^\x - x*|2 - 



i2 1 



(34) 



whereas in region II, under the same conditions, the function u{x) can be expressed as 



Ujj(x) ~ — —{x-x,,) i"^— exp 



'-a^{x - x^Y- 



B2 exp 



-a^x - a;,)2 



(35) 



There is no factor 1/2 in front of the coefficient B2 be- 
cause the expansion of the function Bi differs from the 
expansion of Ai not only in the sign of the argument of 
the exponential but also by a factor 1 /2. Now we evaluate 
in region I, the approximate solution (|29|l for the poten- 
tial ()32|l . The integral in the exponent can be written as 



PX PX^ pX i*X 

/ t^(2/)d?;= / uj{y)dy+ / uj{y)dy = ^ + / ujiy)dy. 

J Xi J Xi J J x^ 

(36) 

The frequency H32|l is used in the second integral only, 
assuming that x is not too far away from x*. The quan- 
tity $ is just a number and its calculation would require 
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the knowledge of 0-1(2;) in the whole region I. However, it 
does not enter the final result and therefore we are not 
interested in its value. We find 

— Aor^\x — exp^— |a; — x^\^ + , 

(37) 

from which we deduce 



be done is to establish the link between C+, C_ and 
82- This can be done by applying the very same 
reasoning in region II. One introduces ^, defined by 



(39) 



Bi = iB2, B2 - A^i 



-l/6pj(<I>-7r/4) 



(38) 



Therefore, we have reached our first goal, i.e. connect 
the two coefficients B\ and B2 to A. What remains to 



and the WKB solution given in Eq. (|31|l reduces to the 
expression 



w„(x) 



'{x — x^) 4 < C+e* exp 



-a'i{x - x,)2 



exp 



-a^ [x - x,)^ 



(40) 



By comparison with Eq. 
C_ are obtained 



3, the expressions for C+ and 



C+ = 



B' 



•2^ 



2 Q,l/6g-*+J7r/2 



B' 



C_ = ^a^/^e*. (41) 



Finally, putting everything together, we arrive at 

Stt 



2*1 



2* 



(42) 



If we only consider the growing mode given by the branch 
proportional to C_, which is a good approximation for 
X ^ X*, one can deduce the power spectra at x = X{. Us- 
ing Eqs. lO, one obtains the dimensionless power spectra 
within the WKB approximation: 



-Pc = 



k 



2^ 
e s 



k 



(l-ei)Kr 
e^'i'T 

(i-£i)KI' 



(43) 



(44) 



where all the quantities in the above expression are eval- 
uated at X = Xf. We know that in the absence of entropy 
perturbations and decaying modes, C is conserved on su- 
perhorizon scales |6|, lll| . Similarly, h is conserved on su- 
perhorizon scales, if the decaying mode is absent. Thus 
we know that both spectra should be time-independent. 
However, within the WKB approximation one cannot see 
this from the above expressions. In the following sections 
we demonstrate this property explicitly for power-law in- 
fiation and slow-roll inflation. 

We see that the calculation of the spectrum for any 
inflaton potential has been reduced to the calculation of 
a quadrature, namely ^ . We are also in a position to 
calculate the spectral indices and their running. The 
spectral indices read 



d^ 

1 = 3 + 2 -'^ 



dlnfc 



= 3 + 2 



d*. 



k=k. 



dlnfc 



(45) 



and the runnings are given by 
. d2*. 



(dlnfc)2 



d2*. 



fc=fc. 



(dlnfc)2 



(46) 



k=k. 



The general expressions for the WKB power spectra, 
spectral indices and their runnings, constitute the main 
result of this work. In order to discuss the power and po- 
tential of the WKB approximation, we now turn to two 
cases where the quadratures and ^'^ can be calculated 
explicitly: power-law and slow-roll inflation. 



IV. APPLICATION TO POWER-LAW 
INFLATION 

We first apply the formalism developed in the previous 
section to power-law inflation. In this model, the scale 
factor is given by 



(47) 



where /3 < —2. This scale factor is a solution of the Ein- 
stein equations in the case where the inflaton potential 
is cx exp((p/mpj). The case /3 = — 2 is special, since 
it corresponds to the de Sitter space-time. The Hubble 
radius is then constant and given by £0. In the general 
case, the horizon flow functions parameters are easily cal- 
culated and read 

2 + ^ - . ,4a, 



ei 



1 + /3' 



e„ = 0, 71 > 1 , 



The exact power spectrum can be found and is given by 



k=k. 



£q Trei 

Vh = $-/(/3)fc2^+^ 



(49) 



(50) 
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where = m ^ is the Planck length and the function 

PI PI o 

/(/3) is given by 



r(-/3-i/2) ' 

2/3+1 



(51) 



In the previous expression, F denotes Euler's integral of 
the second kind. We have /(/3 = —2) = 1. The case 
(3 = —2 is singular, since ei = and the expression of 
the scalar power spectrum blows up. This case should be 
considered separately and one can show that there are 
no density perturbations at all in the de Sitter space- 
time pj| . The spectral indices and runnings can also be 
calculated very easily. One finds — I = =2/3 + 4 
and ttg = a,^ =0. Exact scale-invariance is obtained 
for (3 = —2, which is precisely the singular case evoked 
above. 

Let us now apply the WKB approximation to the case 
under consideration. First, suppose that we want to use 
it before doing the transformation l|24() . In this case, we 
have up {rf) = k"^ — I3{P + 1) I rj^ . In the limit where k-q 
goes to zero, it is easy to see that 



Q 



1 



4|/3(/3+l)| 



(52) 



In general the condition \ Q/uj^\ ^ 1 is not satisfied unless 
—(3 goes to infinity, an unreahstic situation. For /3 ~ —2, 



one recovers that \Q/u 



1/8. This is in agreement 



with the calculation of jQ/w^j performed previously in 
the slow-roll approximation since this approximation is, 
very roughly speaking, an expansion around the scale- 
invariant solution. Equation (|52|l shows that the WKB 
approximation cannot be applied to the usual equation of 
motion without first performing the transformation (|24|l . 
This conclusion is reinforced by the following argument. 
The solution on superhorizon scales is exactly known and 
can be expressed as 



dr 



(53) 



where Ci and C2 are two arbitrary constants. The first 
branch is the growing mode whereas the second one is 
the decaying mode. For power-law inflation this gives a 
mode proportional to /i cx l??]^"'''^ and one proportional 



to /i cx |?7|~'^. One can wonder what the link between 
this solution and the WKB solution is in region II. On 
superhorizon scales, we have ~ [/3(/3 + 1)]^^^/|?7|. As 
a consequence, the mode function /i^^B the form 
A^wKB « |^|i/2±[/3(/3+i)]'/"_ This is not the correct be- 
haviour given by the previous exact solution. The cor- 
rect solution is obtained by replacing the factor I3{(3 -\- 1) 
by (/?+ 1/2)^, i.e. exactly the remark made by Kramers, 
Young and Uhlenbeck (see the Introduction). The rea- 
son is very clear: the effective power-law potential has 
exactly the same shape as the one describing the motion 
in a central-force field, provided that the identification 
t ^ (3 \s made. Let us also remark that the form of 
is the same as the wave function for small r, as ex- 



pected from the previous considerations. This confirms 
that the WKB approximation cannot be used naively on 
superhorizon scales. 

The cure is to use the transformation given in Eq. H24|l . 
After this, the effective frequency takes the form 



ulix) ^ u;l{x) = (1 + /3)V 



-2x 



(54) 



In the superhorizon limit, a; — > cx), one now has Q — s- 0, 
indicating that the WKB solution gives a good approx- 
imation. In addition, the mode function now has the 
correct time dependence. Indeed, since J^' |w|d?/ = 
\(3 -\- l/2|(a;f — x), returning to the dependence in confor- 
mal time gives cx Taking into account the 

link between u and /x, one finally arrives at cx l??!^^'' 



or cx 



i.e. the growing and decaying modes of 
Eq. (|53() in the case of power-law inflation. This means 
that, in this case, the superhorizon solution is nothing 
but the WKB solution. However, this need not be true 
for more general models. A change of variables, which 
allows us to come back from x to the conformal time 77, 
shows that the quantity 5* is given by 



d{kri) 
krj 



(55) 



where rj^ can now be expressed as kij^, =13 + 1/2. This 
integral can be performed exactly, for instance with the 
help of formula (2.275.3) of Ref. (1J|. One obtains 



^' = - 




\(3 + 1/2| - ^-kSf + {f3 + l/2y 



\(3 + 1/2| + + 1/2)2 



(56) 



(57) 



2\f3+l/2\ 

where in the last expression we have considered that \kr][\ <C 1. Inserting this result in Eqs. H43|l and (|44|l . one flnds 



+4 



7'/. = f-5(/3)fc^''+^ m 



2e2/3+i 



(2/3+1)2/3+2 



(58) 
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These expressions are the main result of this section. They should be compared to the exact result given in Eqs. H49|) 
and H5U|I . As expected, they are time- independent. 

I 



The first observation is that the shape of the power spec- 
tra is exactly reproduced by the WKB approximation. In 
particular, this means that the spectral indices and their 
runnings (as well as the runnings of the runnings, and so 
on) are predicted exactly by the WKB approximation. 
The amplitude is given by the function g{f3), whereas 
in the exact case, it is given by the function f{f3); see 
Eq. (inUl- We have g{f3 = -2) = 18 x e^^ ~ 0.896, a 
value that should be compared with /(/3 — —2) = 1. In 
the limit where |/3| goes to infinity we can use the Stir- 
ling formula, T{x) ~ \/27ra;^~"'^/^e~^ to approximate 
and one finds that lim^^_oo /(/?) g{P), i.e. the WKB 
approximation gives the exact result in this limit. It is 
interesting to notice that this limit corresponds to going 
away from scale invariance, i.e. to "fast-roll" inflationary 
models (the kinetic energy of the inflaton is of the same 
order as its potential energy). Therefore, the power-law 
inflation with \(3\ 3> 2 is an explicit example of a model 
where the WKB approximation gives an accurate result 
in a regime where the slow-roll approximation cannot be 
used to compute the power spectrum. 

Another remark is in order. The fact that the WKB 
approximation is better for large values of |/3| than for 
small values of |/3 -I- 2| (scale invariance) might seem to 
be in contradiction with Eqs. (|27|l and H28|l . which in- 
dicate that \Q/uj^\ decreases as the slow-roll parameters 
approach zero (scale invariance). In fact, there is no con- 
tradiction because the WKB criterion only measures the 
accuracy with which the mode function is predicted at a 
given time, without taking into account the error com- 
ing from the matching procedure, while the error in the 
power spectrum is a combination of these two sources of 
error. 

The conclusion is that the WKB amplitude of the spec- 
trum is always a good approximation, see Fig.^ although 
the precision of the prediction for the amplitudes in the 
case of almost scale-invariant inflation is not as good as 
in the slow-roll approximation 0, 0, 0, . Contrary 
to the latter, the spectral indices are predicted exactly 
by the WKB approximation, — 1 = = 2P + A. This 
is a big advantage if data sets of large dynamic ranges 

I 



Vcik) = 



where C = -f^+\n2-2 ~ -0.7296, 7^, ~ 0.5772 being the 
Euler constant. All quantities are evaluated at 77*, which 



are investigated. 



V. APPLICATION TO SLOW-ROLL INFLATION 

Let us now compare the WKB approximation with the 
(leading-order) slow-roll approximation. A comment is 
in order here on the meaning of the phrases "slow-roll 
conditions" and "slow-roll approximation" . If the slow- 
roll conditions hold, all slow-roll parameters are small. 
This does not specify the approximate scheme that one 
should use in order to calculate the power spectra; in 
particular it does not imply that the slow-roll approxi- 



1.00 




0.88 - 

0.85 I I I I 

-14 -11 -8 -5 -2 

FIG. 1: Evolution of the ratio g(l3)/f{l3) with /3. In the limit 
|/3| — > 00, the ratio tends towards 1. This can only be guessed 
from the figure because the convergence rate is slow. 



mation has to be utilized in this situation. In this section, 
we will determine the power spectra computed with the 
help of the WKB approximation in a situation where the 
slow-roll conditions are valid. Having clarified this point, 
let us first quote the result of the slow-roll approxima- 
tion IHifl: 



(59) 
(60) 

I 

is the moment of time when fc* = aH{N^). The scale fc* 
is called the pivot scale. It has been shown in Ref. Q 



nm'i. 



1 - 2(C + l)ei - Cea - (2ei + £2) In 



1 - 2(C+ l)ei - 2eiln 
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that the corresponding time ry, has a nice interpretation: 
it is nothing but the time at which the matching between 
the region where aH ~ k and the superhorizon region is 
performed. It is important to reahze that, according to 
the previous definition, 77* is not a function of k, i.e. is 
the same for ah scales. Changing the time of matching 
is therefore equivalent to changing the pivot scale. This 
means that the pivot scale must be in the range of scales 
that are of astrophysical interest today but otherwise can 
be chosen freely. It has been demonstrated in Ref. 0| 
that the location of the pivot scale can be optimized to 
reduce the slow-roll error in the power spectrum. 

From Eqs. H59I) and H60() . we notice that close to the 
de Sitter inflation the slow-roll approximation predicts 
the amplitude for both density and tensor perturbations 
very well. The accuracy of the slow-roll approximation 
improves as one approaches scale-invariance, since this 
corresponds to e„ — > (for power-law inflation this cor- 
responds to /3 — > —2). A priori^ it gives an exact result 
for exact scale invariance but then the linear approxima- 
tion blows up, as can be seen from Eq. (|59|l . when e — > 0. 

In contrast, the WKB approximation does not give an 
exact result for the amplitude in the scale-invariant limit 
(see the discussion of the limit (3^—2 above), but it 
is able to give a reasonable approximation in the regime 
where slow-roll is violated (ei 1). The spectral indices 
predicted by the slow-roll approximation at first order 
are 



frequencies read 



1 



-2ei - 62 



(61) 



These estimates are not exact, contrary to the WKB re- 
sult in the case of power-law inflation, but provide good 
estimates close to the scale-invariant situation. 

Let us now evaluate the WKB power spectra for the sit- 
uation of slow-roll inflation. For that purpose we expand 
(|43|l and H44|) at flrst order in the horizon flow functions 
and compare the result with (|59|l and H60() . The effective 



(62) 
(63) 



and the WKB power spectra become 



p,. = M!(-M'|(i-|.i)«»', 

TTTO^ 3 V 3 



(65) 



where we have made use of the relation 1/[(1 — ei)aH] = 
— ry[l -I- 0(e^)]. The evaluation of the integrals in the 
argument of the exponential functions gives 

^{m.V*;k) = ( ' W{y;k)\Ay (66) 



with 



^ k^T^+b{e,5,0—. (67) 

4 T 



6g(e„) = 3ei-f-e2, Ki^n) = ^^i , (68) 



r?*(5) 



1 /9 



fc V 4 



+ 6(en). 



(69) 



At flrst order, we may consider the horizon flow func- 
tions as being constants. This is due to the fact that the 
derivatives of the horizon fiow functions are quadratic in 
these parameters. The integration can be performed as 
for power-law infiation in the previous section. The final 
result is 



Vcik) - _^(18e-3) [1 _ 2{D + l)ei - De^ - (2ei + £3) ln(-fcr?f) + ©(e^ )] , (70) 

^ PI 

Vh{k) = i^(18e-3) [1 _ 2{D + I)ei - 2ei ln(-fcr;t) + 0{el)\ , (71) 
PI 

with D = ^ — ln3 « —0.765, a reasonable approximation of C from above. All the quantities in the above equation are 
evaluated at tj = rjf. In order to compare this result with the slow-roll expressions H59|l and H60|) we have to establish 
a link between H{Nf) and H{N^,), ei{Nf) and ei(iV*), and so on. At the level of a first-order expansion in the horizon 
flow functions, we can consider ei and €2 to be constant everywhere, except in the factor 1/ei in (|59(l . Performing a 
linear expansion of ei around iV* we find 

^ ^ -[1 - e2{N^)AN + . . .]. (72) 



ei(iVf) ei(iV,)' 
In a similar manner we find 

H^{Ni) ^ H^{N^)[1 - 2ei{N^)AN + . . .]. (73) 
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We eliminate —rji in the logarithm by expressing it in terms of l/{aH){Nf). Corrections containing ei inside the 
logarithm can be omitted since they contribute at higher order only. Finally, the logarithm can be written 



ln(— /c77f) — In 



k 



{aH){N, 



In 



iaH)iNi) 



= In ( A ) _ A7V 



(74) 



where now the dots denote terms of order ei and we have used the pivot scale fc» = {aH){N^,). Putting everything 
together we find that all AiV cancel, and the result simply reads 



Vcik) 



neim: 
16iJ2 
■nm? 



-(18e"3) 



1-2{D + l)ei - De2 - (2ei + £3) ln( ^ ) + 0{el) 



(18e-3) 



1 - 2(D + l)ei - 2eiln 



0{el) 



(75) 
(76) 



which can now be compared with the prediction from the 
slow- roll power spectra (I59|l and Ht)0|) . 

It is seen immediately that the spectral indices coincide 
with the slow-roll result (|61|l . As already expected from 
the discussion of power-law inflation above, the spectral 
indices are predicted by the WKB approximation very 
well, the amplitudes only at an accuracy of 10% (coming 
from the factor 18e~'^). This is certainly not good enough 
for future experiments, but it is well known that WKB 
approximations can be systematically improved. 

VI. CONCLUSIONS 

A new method to predict inflationary power spectra 
has been presented, based on the WKB approximation. 
One of the main advantages of this new scheme of approx- 
imation is that it can give a good approximation even if 
the slow-roll conditions are violated, whereas, up to now, 
all of the methods in the literature were subject to this 
limitation We have checked this fact explicitly on 

the example of power-law inflation. In addition, we have 
also tested the method for the case of slow-roll inflation 
and find errors below 10% for all quantities considered. 
The prediction of the spectral indices is exact in the case 
of the power-law inflation. Assuming that the conclu- 
sions obtained from the power-law case are also valid in 
a more general context, the WKB approximation appears 
as an efficient method, especially for predicting the spec- 
tral indices. The main restriction comes from the ampli- 
tude, which is not predicted with a similar accuracy as 
the slow-roll approximation close to scale invariance. 

A possible disadvantage of the method is that it is non- 
local in the sense that an integration is involved from 
the time of horizon crossing until some time when all the 
modes of interest are on superhorizon scales. This means, 
in contrast to the slow-roll approximation, where a cer- 
tain behaviour of the background is assumed, that one 
needs to know how the Hubble flow functions evolve with 
time in principle. In the example of slow-roll inflation it 
has been shown that one essential step in the WKB ap- 



proximation is to link the predictions at some time rj^ to 
the parameters relevant at the time of horizon crossing. 

Many improvements of the results presented in this 
article are possible. As already mentioned above, the 
WKB expansion can be pushed further, e.g. by going 
beyond the linear approximation for the joining of the 
super- and subhorizon WKB solutions. This will reduce 
the error in the amplitude, therefore avoiding the main 
shortcoming of the new method. Another point is that 
the standard slow-roll approximation (Bessel approxima- 
tion) is valid only when the slow-roll parameters are con- 
stant. Therefore, it cannot be pushed beyond the first 
order. Recently, Stewart and Gong ^T§\ have proposed 
a method, based on Green fimctions perturbative cal- 
culations, which is free from this problem. They have 
calculated the inflationary power spectra up to the sec- 
ond order. The WKB approximation is also free from the 
above limitation and therefore allows a determination of 
the spectra beyond the leading order. All these questions 
are currently under investigation 17]. 

Note added: While the present article was being com- 
pleted, a new paper |l8j | appeared where the inflationary 
power spectra are also calculated using a new method of 
approximation, namely the "uniform approximation" ac- 
cording to the terminology used in that article. Some of 
the equations obtained in bear a close resemblance 
to the results obtained in the present work as revealed 
by a comparison of Eqs. and with Eq. (14) of 
Ref. After having released the first version of this 

paper, we became aware of an earlier attempt to use the 
WKB approximation for cosmological perturbations by 
Nagasawa and Yokoyama ■ 
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